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Abstract 

We describe a numerical simulation of the evolution of an S '3 cosmic string network 
which takes fully into account the non-commutative nature of the cosmic string fluxes 
and the topological obstructions which hinder strings from moving past each other 
or intercommuting. The influence of initial conditions, string tensions, and other 
parameters on the network’s evolution is explored. Contrary to some previous sug¬ 
gestions, we hnd no strong evidence of the “freezing” required for a string-dominated 
cosmological scenario. Instead, the results in a broad range of regimes are consistent 
with the familiar scaling law; i.e., a constant number of strings per horizon volume. 
The size of this number, however, can vary quite a bit, as can other overall features. 
There is a surprisingly strong dependence on the statistical properties of the initial 
conditions. We also observe a rich variety of interesting new structures, such as light 
string webs stretched between heavier strings, which are not seen in Abelian networks. 
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1 Introduction 


A generic feature of many spontaneously broken gauge theories is the existence 
of topological solitons, such as strings (or flux tubes), domain walls and monopoles. 
Many grand unihed models predict the formation of such defects during a cosmolog¬ 
ical phase transition!^' The defects could be interesting as a potentially observable 
signature of the symmetry-breaking pattern, and could also have important conse¬ 
quences for the evolution of the universe. Domain walls and monopoles, if they are 
stable and occur without cosmic strings, generally are not considered phenomeno¬ 
logically viable. Domain walls tend to produce density perturbations that are much 
too strong, whereas monopoles, if they exist at all, are predicted to form in great 
abundances which are incompatible with current observations. Among the possi¬ 
ble categories of stable topological defects, string-like defects, or cosmic strings, are 
considered the least disastrous for cosmological models, and may be useful from a 
model-building point of view!^' Among the potential applications, it has been pro¬ 
posed that the gravitational effects of either inhnite cosmic strings or closed loops 
of string may serve as sources of density perturbations leading to galaxy or cluster 

r .• [5.61 

formation. 

A persistent string network could conceivably also have profound effects on the 
evolution of the universe due to its bulk energy density, quite apart from the effects 
of fluctuations. If a network of strings becomes frozen so that strings are hxed in 
comoving coordinates, then they will be stretched by the expansion of the universe. 
If the network is thought of as composed of a fixed number of segments, the number 
of segments per unit volume will be proportional to a~^ while the length of each 
segment grows as a due to stretching, and so the total energy density will scale as 

a ^ X a a 

where a is the scale parameter representing the size of the universe. The energy 
density of non-relativistic matter, on the other hand, scales as a~^, and as the universe 
expands the energy in strings will grow relative to that of matter until it eventually 
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dominates. As explained below, such a frozen network is not a typical outcome for 
the types of strings that have been studied to date. Instead, there are energy loss 
mechanisms that allow strings to be progressively destroyed, while those that survive 
continue to move relativistically, their total energy scaling in the same way as that 
of matter. It has been suggested, however, that non-Abelian strings might behave 
differently from Abelian ones, and might indeed lead to a string-dominated universe. 
A universe dominated by very heavy strings is not likely to be a viable model, but a 
cosmological model with a fairly recent transition to a string-dominated phase with 
comparatively light strings has some desirable properties. Strings could serve as 
an interesting form of dark matter: giving density parameter hi ~ 1 as required by 
inflation, while ^matter < 1, and possibly resolving the apparent discrepancy between 
estimates of the age of the universe from the expansion rate and from stellar ages. 
By modifying the equation of state of the universe, they could mimic some of the 
effects of a cosmological constant. There has been a recent revival of interest in such 
a scenario, and testing its consistency is one of the chief motivations for the work 
described in this paper. 

The evolution of Abelian cosmic strings has been studied extensively, and we 
review some of the salient points here in order to point out contrasts with the non- 
Abelian case. The simplest and best understood type of cosmic strings are those 
which occur in the Abelian Higgs model and are classihed by their integer winding 
number (an element of Z). We shall refer to these as Z strings.* Assuming that 
only the strings with unit winding number are stable (type II strings), they can be 
formed either as inhnite strings or closed loops. Monte Carlo simulationsindi¬ 
cate that the inhnite strings constitute a majority (63-80%) of the total length of 
the network initially formed by a phase transition, with the remainder comprising a 
“scale-invariant” distribution of closed loops. A very important dynamical process 
of these strings is intercommutationj^’^^’^^' the process in which two colliding strings 

* These are often referred to in the literature as U{1) strings since they result from the complete 
breaking of a C/(l) gauge group, but it seems more logical to call them by their topological 
classification, as is usual for other types of strings. 
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reconnect as in Figure 1. Two infinite strings intercommuting with each other twice 
(or a single string intercommuting with itself) can form new closed loops from pieces 
of the inhnite strings, and this is thought to be the principal mechanism for the 
destruction of inhnite strings. 

Dimensional analysis and energy conservation arguments provide a plausible 
picture of the strings’ evolution: Inhnite strings are destroyed by intercommutation 
and thus lose energy to the network of loops. Loops, in turn, shrink as they lose 
their energy to gravitational radiation, and they may also split into smaller loops by 
intercommutation. Any loop will eventually contract to a point and be destroyed. 
The system reaches a dynamical steady state, or self-similar evolution, characterized 
by a single evolving length scale which is the size of the cosmological horizon. The 
rate at which new strings appear within the growing horizon is balanced by the rate 
at which they are destroyed. An approximately constant number (order unity) of long 
strings stretch across a Hubble volume at a given time, and intercommutations result 
in the formation of a similar number of new closed loops per Hubble volume, which are 
destroyed within a Hubble time. While there is disagreement over details, numerical 
simulations have supported this general picture. All results are consistent with a 
scaling behavior in which the energy density in strings remains a small, hxed fraction 
of the matter density and there is no string domination. 

We emphasize that intercommutation is crucial to this picture of the strings’ 
evolution. Since non-Abelian strings in general cannot intercommute, we might expect 
different behavior for such strings. 

Somewhat less attention has been devoted to the evolution of branched networks, 
in which several strings may join at a vertex. Branched networks occur when a U{1) 
gauge group is spontaneously broken to with N > 3, or when the unbroken 
group is a non-Abelian discrete group. Among the work that has been done is that 
of Vachaspati and Vilenkin who considered a network of strings, which have 
the novel feature that three strings may intersect in a vertex. Z 3 strings tend to 
form an inhnite network of vertices connected by string segments, with very few 
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Figure 1: Intercommutation of colliding strings. 
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closed loops. It had previously been speculated that the nodes in a branched network 
could settle to equilibrium positions, thus causing the network to freeze as a string- 
dominated universe requires, but the simulations of [14] indicated otherwise. Instead 
of reaching an equilibrium, the nodes pull together and annihilate, steadily reducing 
the number of nodes and strings. The annihilation of vertices leads to a self-similar 
scaling behavior, as long as the nodes are able to come close enough to each other to 
annihilate. The phenomenological consequences of their model are rather similar to 
those of Z strings: The number of vertices and string segments per horizon volume 
remains roughly constant, and the energy density of the network is a small constant 
fraction of the matter density. The strings never relax to an equilibrium, but continue 
to move with relativistic transverse velocities, following a self-similar evolution pattern 
much like that of Z strings. 

The purpose of this paper is to explore the consequences of a network of non- 
Abelian strings. Such strings are known to exhibit a number of exotic types of in¬ 
teractionsParticularly signihcant is the fact that when two non-Abelian strings 
cross each other, they cannot generally intercommute, nor can they pass through 
one another without forming new vertices and becoming joined by a new segment of 
string. Linked loops of string cannot usually become unlinked, and vice versa. One 
might expect that this would inhibit the decay of a cosmic network by obstructing the 
removal of string segments. If new strings are continually being formed through string 
collisions, their energy must come from the already existing strings, and an equilib¬ 
rium with the strings’ transverse velocity damped out might seem a more likely hnal 
state. Some evidence for slowing down of the network’s destruction was reported very 
recently by Pen and Spergel in a class of models with non-Abelian global strings. 

In this paper, we describe a numerical simulation of a network of Ss strings. 
Unlike the authors of reference [8], we consider gauge strings, which have no mutual 
long-range interactions. Our interest is in understanding the qualitative nature of 
the fate of a non-Abelian string network. Do collisions result in a non-diminishing 
or rapidly increasing number of strings? A static equilibrium state which is expected 
to be conformally stretched with the universe’s expansion? Does the network instead 
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decay rapidly into finite networks and closed loops? Or does a dynamical self-similar 
evolution emerge, as in [14]? Another question which we hope to illuminate is: which 
processes play the most crucial roles in the network’s evolution? The importance 
of intercommutation to the evolution of Z strings led to fundamental study of the 
dynamics of intercommutation; likewise it is hoped that the results obtained here will 
suggest which aspects of non-Abelian string dynamics are ripe for closer examination. 

Our method of simulation is directly inspired by that of reference [14]: we generate 
initial conditions from a lattice Monte Carlo simulation and then evolve the network 
according to a highly simplified model of string dynamics which we hope captures 
the essential features of a string network losing energy. We find hints of some quite 
interesting physics in the interplay between the two types of string in our model, and 
a rather surprising dependence of the network’s behavior on the initial conditions. 
Concerning the string-dominated cosmological scenario, we reach somewhat different 
conclusions from [ 8 ]. For a wide range of conditions, the network’s density follows 
power laws much like the ones arising in Abelian networks. Non-abelian effects can, 
indeed, slow down the network’s decay in the sense that they change the coefficients of 
the power laws, but it appears as if the slowing down to a stable equilibrium happens 
only under special circumstances, if at all. 

The remainder of this paper is organized as follows: Section 2 provides a brief 
summary of the properties of non-Abelian strings which bear on this simulation. We 
discuss some of the subtleties inherent in the description of non-commuting magnetic 
fluxes, and the necessity of a (gauge fixing) convention to resolve these ambiguities 
and allow the comparison of fluxes of strings. Most importantly, we explain why 
two colliding non-Abelian strings cannot, in general, cross or intercommute without 
forming a new segment of string. Section 3 describes the particular S '3 model which 
we have chosen to simulate. S' 3 , the permutation group on 3 elements, was chosen 
as the gauge group because it is a simple non-Abelian group which exhibits all of 
the important general characteristics of non-Abelian strings. Another motivation for 
this choice is that S '3 contains Z 3 as a subgroup, allowing instructive comparisons 
with the Z 3 network simulations of ref. [14] . Section 3 also describes our procedure 
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for simulating the network’s dynamical evolution, giving enough details to allow an 
understanding of the results. Additional technicalities of the procedure are relegated 
to the appendix. Section 4 describes our procedures for generating initial conditions, 
and summarizes features of the networks these procedures generate. We use two 
different Monte Carlo algorithms which generate initial networks with somewhat dif¬ 
ferent statistical properties. The network’s evolution turns out to have a surprisingly 
strong dependence on the initial conditions. Section 5 presents results of the dynam¬ 
ical evolution simulation, exploring the influence of a number of different variables 
including the initial conditions of the network and the ratio of string tensions. These 
results are compared with those for a Z 3 network, which is Abelian. In section 6 , we 
present our conclusions and suggest directions for future work. 

An appendix describes our procedure for keeping track of string fluxes during 
the simulation and for establishing them from the lattice Monte Carlo procedure, 
covering details not included in section 3. The implementation of non-Abelian fluxes 
in a simulation presents a rather difficult problem in its own right. A careful gauge- 
hxing procedure is required, and some of the subtleties that arise are of interest from 
a held-theoretic point of view. The algorithm has been described in greater detail in 
reference [18]. 

In this discussion, the strings will be considered as classical objects with well- 
dehned fluxes (after a gauge has been fixed). We will not consider quantum- 
mechanical effects such as Cheshire charge.'^*^^ 

2 Vortices and strings in a non-Abelian discrete gauge theory 

In this section, we review very briefly the dehnition of non-Abelian vortices and 
strings and some of the properties which are important for the current simulation. 
The formalism used here was developed for vortices in [20] and applied to strings in 

|16]. 

Generically, topological defects of codimension 1 (vortices in two space dimen¬ 
sions, strings in three) occur when a gauge symmetry group G is spontaneously bro- 



ken to a subgroup H such that there are non-contractible closed loops in the vacuum 
manifold G/H. This is formally expressed by saying that t:i{G/H) ^ 1, where tti 
stands for the hrst homotopy group!^’^^^ Each element of tti represents a closed path 
(or, more precisely, a class of closed paths.) Paths can be composed by tracing hrst 
one path and then another, giving the group structure of tti. For example, if a U{1) 
gauge group is broken completely, then H = 1, and ni{G/H) = Z, the group of the 
integers. Closed paths in the vacuum manifold are classihed by an integer winding 
number. If G is simply connected, then 7ii{G/H) ~ no{H), where ttq is the set of dis¬ 
connected components, so that strings occur whenever H is discrete or has components 
disconnected from the identity. We will subsequently assume G is simply connected 
unless otherwise stated. A further simplihcation is when H is discrete, in which case 
no{H) H. If 71q{H) is non-Abelian, then the composition of paths depends on the 
order. Hence the fluxes of strings will be non-commuting group elements, and that is 
what is meant by non-Abelian strings. 

The breaking of an underlying G gauge theory to a discrete group H leaves no light 
propagating gauge helds: At low energies, in any simply connected region without 
defects, the gauge held A“ is pure gauge. However, when string defects are present, 
the region of true vacuum is not simply connected: iiisTZ = M. — {T*}, where {D} is 
the union of all defect cores (regions of non-vacuum) and M. is the spatial manifold 
on which the defects exist. Each string gives rise to a class of noncontractible closed 
paths in AI — {D} which encircle the string. The hux enclosed by any closed loop T 
is a group element dehned as a path-ordered exponential of the gauge held: 


fiux = Pexp{j) A-di). (1) 

r 

For any T within TZ this must be an element of H. This is because the Higgs held 
is covariantly constant throughout TZ and so the transformation that results from 
parallel transport around a loop must leave the Higgs held invariant. 

The hux of a cosmic string is dehned by the above exponential along a path 
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which winds around the string. In a non-Abelian theory, this definition of the flux 
is not gauge invariant, and may depend on the point at which the path begins and 
ends. However, the flux through any loop which does not enclose a string is necessarily 
trivial. A corollary of this fact is that two closed loops which share the same beginning 
and ending point xq, and can be continuously deformed into each other, have the 
same flux. The relevant structure for the description of the system of defects is the 
fundamental group or first homotopy group — {T)},a;o), defined with respect 

to a basepoint xq. Each string is associated with a generator of the fundamental 
group. Once xq has been (arbitrarily) chosen, the fluxes of all closed paths (and of all 
strings) are specified by a homomorphism from — {T*}, xq) into H ~ 7ri(G/ff). 

The only remaining gauge freedom is a global one. However, there is a considerable 
amount of ambiguity in what we mean by “the flux” of one particular string: an 
arbitrariness in how exactly the set of generators is chosen for the homotopy group 
7 ri(Al — {Zi)},a;o). In figure 2, for example, there are two loops, both beginning and 
ending at xq, both enclosing the same string without enclosing any others, which are 
nonetheless representatives of different homotopy classes (and consequently may be 
associated with different fluxes). An intervening string prevents one path from being 
continuously deformed to the other. The fluxes associated with the two different paths 
may differ through conjugation by the flux of the other string. In an Abelian theory, 
conjugation is trivial; not so in a non-Abelian one. It follows that fluxes cannot 
meaningfully be compared (say, to determine if they are the same) if the paths used 
to define those fluxes pass on opposite sides of some other string. Comparisons must 
be made using “nearby” paths. 

Since the flux of the same string may described, depending on convention, by 
distinct conjugate elements of H, it follows that strings whose fluxes are in the same 
conjugacy class must be degenerate in tension. It is not true, however, that all fluxes 
in the same conjugacy class are the same. If a and b belong to the same conjugacy 
class, it is by no means guaranteed that and ab also are conjugate, or that ab = ba. 
Distinctions among elements within the same conjugacy class can have important 
consequences in any situation where the product or commutator of two fluxes is 
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relevant, as is the case when two strings collide. 

The path dependence of the flux of a string implies an important fact: two strings 
with non-commuting fluxes cannot pass through each other without forming a new 
segment of string whose flux is the commutator of the fluxes of the two original 
strings. Penetration without the formation of a new string would violate flux con- 
servationl^^^This is illustrated in hgure 3. Non-commuting strings also cannot inter¬ 
commute. We will be especially interested in the consequences of this entanglement 
process for the evolution of a string network: it might impede the collapse of the 
network. 

3 Our Model and its Dynamics 

S '3 strings 

We consider here a model with unbroken gauge group H = S' 3 , the permutation 
group on three objects. The spectrum of this model will include strings whose fluxes 
are elements of S' 3 . S '3 has six elements in all. The identity e corresponds to the 
trivial permutation. There are three odd permutations (two-cycles or transpositions) 
each leaving one of the three elements invariant and interchanging the other two. 
We may denote these, for convenience, by: ti = {(123) —(132)},t 2 = {(123) — 
(321)},t 3 = {(123) —(213)}. In this notation, ti is the two-cycle which leaves 
the i-th element in the same position. The two non-trivial even permutations are 
the three-cycles, or cyclic permutations, which we denote here by s_|_ = {(123) — 
(312)}, s_ = {(123) —(231)}. In the more conventional cycle notation, we have 
ti = (23), f 2 = (13), t 3 = ( 12 ), s+ = (123), = (132). 

The 2-cycles form one of the two non-trivial conjugacy classes, and the 3-cycles 
form another. Thus our model supports two types of strings, which we shall refer to 
as odd and even strings, or alternately as t and s strings. The three-cycles generate 
a Z 3 subgroup, so that our model contains the Z 3 model as a subset. Three even 
strings may meet at a vertex, just as in the Z 3 model. Another type of junction is 
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Figure 2: The paths a and a' both enclose the same string and no other strings, 
but they cannot be continuously deformed into each other without crossing another 
string. Thus, they represent different elements of the fundamental group 7ri(Af — 
{D},a;o), and so the fluxes associated with them may be different. Specifically, the 
homotopy classes of a and a' are related through conjugation by another generator: 
a' ~ (We follow the usual convention of composing paths from right to left: 

(3a(3~^ means the path formed by traversing first the reverse of jS, then a, then (3. The 
relation ~ represents homotopy equivalence.) The associated fluxes are analogously 
related: a nontrivial relation if the fluxes don’t commute. 
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Figure 3: Attempt to pass two strings through each other. In (A) the flux of 
one string may be defined by either of the paths a or 7 , and that of the other string 
by {3. Let the fluxes associated with a, f3 and 7 be a, 6 , and c respectively. In this 
case, c = bab~^. In general, c ^ a. Now, if we attempt to pass the strings through 
each other, no strings need cross the paths a and 7 , so the associated fluxes will not 
change. But if the strings were to pass through each other freely, as in (B), a and 
7 would be continuously deformable into each other. This is impossible if they have 
different fluxes. In order to conserve flux, the string must branch somewhere and 
be connected to the other by a new string whose flux as defined by path 6 in (C) is 
ca~^ = bab~^a~^. 
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one where two 2-cycle (or odd) strings merge to form a 3-cycle (even) string. Figure 
4 shows the two types of junctions in our model. 

Since each two-cycle is equal to its inverse, oppositely oriented odd strings are 
topologically equivalent. An even string, on the other hand, possesses a natural 
orientation; The flux through a path encircling it with one orientation is while 
it is S- for the opposite orientation. In subsequent hgures, even strings will often 
be denoted by oriented lines, with the string carrying flux s_|_ in the direction of the 
arrow, while odd strings have no arrow, reflecting their lack of orientation. 

Note that the parity is conserved at any junction: i.e., an odd string entering a 
junction must be matched by another odd string leaving. In this sense, odd strings 
can never end, even if they change their flux at a junction. Even strings, on the other 
hand, may end at a junction. It is helpful to view the network as being composed 
of two interacting subsystems. One subsystem consists of inhnite or closed t-strings 
(like Z strings, they have no free ends). The other comprises a branched Abelian web 
of s-strings, some of which end on t-strings. The simulation results presented later in 
this paper exhibit an interesting interplay between these two subsystems. 

String Dynamics 

The system we simulate consists of three-way junctions, or nodes, joined by s 
and t strings. The strings are approximated as straight segments between junctions. 
In effect, we are averaging over transverse oscillations of the strings. However, it is 
possible for a string segment to be interrupted by a pair of nodes doubly linked to 
each other as in hgure 5. Doubly linked nodes tend to annihilate fairly quickly, and 
such conhgurations play the role of transient kinks on otherwise straight segments. 

As in [14], the nodes are assumed to undergo damped motion under the influence 
of string tensions. The energy loss leads to the shortening of strings. As an approxi¬ 
mation for the energy loss of a real string network, the model of damped motion of the 
nodes is most realistic if one supposes that the string junctions are monopoles which 
carry some unconhned magnetic flux, as in a model with the symmetry-breaking 
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(B) 


Figure 4: String junctions in the S '3 model. (A) Two possible sss junctions: 
three strings with the same flux, s_|_ or s_, emanate from the node. (Or two s_|_ 
strings merge into a single s_, etc.) (B) One of the class of stt junctions: Two t- 
strings merge into an s-string. Fluxes are defined with respect to a basepoint xq 
by the paths shown. Here, as in many subsequent figures, an s-string is drawn as 
an oriented line. The string carries flux s+ in the direction of the arrow: i.e., a 
counterclockwise path around the arrow encloses flux s_|_ 
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Figure 5: Doubly linked nodes. A pair like this can interrupt an otherwise 
straight string or act as a kink. 
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pattern 


S 3 X (SU(3) X SU(2) X i7(l)Bj,) 

- D -’ 

where D is a discrete factor divided out so that the monopoles at string junctions 
may carry electromagnetic U{1) charge. Their magnetic charges should then result 
in radiation damping.* 

Our simulation proceeds in discrete time steps. During each time step, each node 
is moved by a displacement proportional to the vector sum of all tensions acting on 
it. This type of evolution corresponds to damped motion: Force oc Velocity. The 
constant of proportionality is a parameter which may be absorbed into the size of the 
time step. Thus, in appropriate units. 

Ax = At E Trllr, ( 2 ) 


where At is the time step, n,. is the unit vector along the direction of the r-th string 
connected to the node, and Tr is the magnitude of that string’s tension, r runs from 
1 to 3 for the three strings that meet at each node. Since our model has two types 
of string with possibly different tensions, the ratio of these two tensions is a variable 
parameter of the simulation. When we present the results in section 6 , we use units 
such that the lattice spacing of the initial-condition Monte Carlo algorithm is 1 and 
the tensions are of order unity. (More specihcally, the f-string tension is normalized 
to 1 while the other tension is varied.) 

Each node is moved in turn. During the motion of a node, the strings attached 
to it may collide with other strings. The annihilation of two nodes is allowed if they 

* A model has been constructed in which topological strings become attached to monopoles 
which also carry other charges. Constructing a model with S 3 strings joined at monopoles 
might be slightly harder, but it is not our main concern here. For a hint of how such a model 
could arise, consider the monopoles that form when an SU (5) group is broken in the familiar 
way to SU{3) x SU{2) x U{1)/Zq. This transition is known to yield stable monopoles with 
SU{3),SU{2), and 17(1) We could imagine a second symmetry-breaking stage in 

which the SU{ 3 )/Z 3 factor is broken down to S 3 in such a way that the resulting srings also 
carry nontrivial flux in the Z 2 center of SU (2). Whenever three such strings join, the resulting 
net Z 2 flux can unwind through a monopole, which has both SU(2) and 17(1) flux. 
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approach each other more closely than a certain distance dann- The procedures for 
handling collisions of strings or nodes are as follows: ^ 

i) Intercommutation: If, in the process of moving a node from its initial to hnal 
position, one of its string segments intersects some other segment, then the fluxes of 
those two segments are compared at the point where the crossing occurs. If the fluxes 
commute, then the two segments may either pass through each other unaltered, or 
intercommute. The probabilities of these two outcomes may in principle be taken as 
an adjustable parameter of the simulation, but we have chosen to let the probability 
of intercommutation be 1 in all cases. It is widely believed that intercommutation 
is generically the more common outcome whenever two cosmic strings cross. In the 
self-dual limit, it can be shown that Nielsen-Oleson strings always intercommute. 
Therefore the choice to set the intercommutation probability to 1 seems a natural 
one. 

Intercommutation may occur in two possible situations: either both strings are 
3-cycle (s) strings, or both are 2-cycles. In the latter case, the fluxes of the two strings 
must in fact be equal. In an intercommutation, the string ends are rearranged in such 
a way as to conserve flux. In the case of two s-strings, there is always only one way 
to rearrange the ends, as shown in figure 6a. A string end carrying flux to the 
point of intersection may not be joined to one carrying the inverse flux s_. When two 
t-strings intercommute, however, there are two possible rearrangements of the ends, 
owing to the fact that a 2-cycle is equal to its own inverse and 2-cycle or t-strings 
consequently have no preferred orientation^^ (£g. 6b.) In the absence of a reason to 
prefer one of these rearrangements over the other, the choice is made randomly. 

When an intercommutation occurs, we join the segments with a pair of “kinks” 

t Some of these procedures differ in minor details from those described in [18]. 
ft Strictly speaking, we can only say that there is no topological reason for a t-string to have a 
preferred orientation. It is possible that the field equations could have two distinct solutions, 
corresponding to differently oriented strings, which are topologically equivalent but can be 
deformed into one another only by surmounting a finite energy barrier. A situation of this 
sort occurs in the global vortices of nematic liquid crystals. This was pointed out to me by J. 
Preskill. 
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(doubly linked nodes as in figure 5) which will later annihilate and allow the segments 
to straighten. By delaying the straightening of the rejoined segments, we prevent the 
highly non-causal cascades of intercommutations which might otherwise occasionally 
occur within a single time step and lead to computational infinite loops and other 
unpleasant consequences. 

ii) Non-commuting collisions (NCC): If two non-commuting strings intersect, 
then it is assumed that they form a new pair of nodes and thus become linked by 
a new segment. This may happen in two possible ways. The two strings may pass 
through each other and become linked by a new segment which stretches between 
them, as shown in figure 7a. We refer to this outcome as the “bridge” configuration. 
The flux of the intervening string segment is uniquely determined by the requirement 
of flux conservation. (The intervening flux of the bridge must always be a 3-cycle, 
as the commutator subgroup of S '3 is Z 3 .) Another possible outcome, consistent 
with the topology, is that the two colliding segments stick together in the “zipper” 
configuration shown in figure 7b. The possibility of collisions resulting in zipper 
formation have recently been mentioned in the context of type I Abelian strings 
Which of these two outcomes is more likely is a dynamical and kinematic issue to 
be addressed in future work. It may depend on the relative velocities and angles 
of the colliding strings. In the present simulation, we consider both possibilities 
separately and examine their consequences: some runs were performed with only 
“bridge” NCC’s, and some with only “zipper” formation. In fact, we find that the 
choice makes little difference for the network’s evolution. 

Whenever two new nodes are created by a NCC, we separate them immediately 
by some small initial distance comparable to the time-step and subesquently allow 
them to move normally under the influence of string tensions. 

iii) Annihilation: When two nodes approach each other within a distance dann 
which is a parameter of the simulation, they are allowed to annihilate. (In this work, 
dann = 0.08 was chosen.) The segment(s) which join the two nodes is eliminated, 
and the other segments emanating from the two annihilating nodes are joined to each 
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Figure 6A: When two three-cycle or s-strings intercommute, there is a unique 
rearrangement which is compatible with the orientations of the strings. 
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Figure 6B: Since t-strings (2-cycle strings) have no orientation, an intercommu¬ 
tation can result in either of two possible rejoinings of the ends. 
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Figure 7: The intersection of two strings whose fluxes do not commute causes 
them to become linked by a new segment in the “bridge” conhguration (A). Alterna¬ 
tively, they may coalesce along part of their length, forming a “zipper” (B). 
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(A) "Bridge" (B) "Zipper' 



Figure 8 : A pair of nodes linked by two strings may annihilate, leaving a single 
string. 
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Figure 9: Annihilation of two sss nodes joined by a single string. There are two 
possible ways to reconnect the strings consistent with their orientation. 
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Figure 10: The two nodes shown here cannot annihilate, because there is no 
consistent way to reconnect the string ends. 
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other. 


Two nodes are able to annihilate only if there is a consistent way to rearrange 
the free string ends (i.e., each string is able to hnd a partner with the same flux). 
Annihilation is always possible if the two nodes are doubly linked as shown in hgure 
8. In the case of double link annihilation, the segment is straightened (or straight¬ 
ened until an obstruction is encountered, such as a collision with another string.) 
Annihilation is also always possible if both junctions are of the sss type, even if they 
are only singly linked. In this case, there are two possible rearrangements of the free 
string ends (hgure 9). One of these two is chosen at random. When two stt-type 
junctions approach each other, on the other hand, there can be at most one consistent 
rearrangement of the free ends, and it may not be possible for the nodes to annihilate 
at all. Annihilation requires that each of the two segments on one side be matched 
with one on the other side carrying the same flux. Figure 10 shows an example of a 
pair of nodes which cannot annihilate because there is no consistent rejoining of the 
string ends. 

If two stt-type nodes do annihilate, it is easily seen that there can never be more 
than one consistent rearrangement of string ends. If two of the outgoing ends are 
s strings and two are t strings, then there cannot be more than one rearrangement 
because each string can only be joined with one in the same conjugacy class. If all 
outgoing strings are of odd type, then all four cannot have the same flux- if they did, 
then the total flux of any pair would be trivial and they would not be connected by a 
segment. Nor may any three have the same flux. It follows that, at best, each string 
end may reconnect with a unique partner. 

When a singly linked pair of nodes annihilates, the segments are hrst rejoined in 
a kink conhguration, which may straighten later. 

Another type of annihilation process, which is the inverse process of bridge for¬ 
mation, is not included explicitly in the simulation but may occur through through a 
multi-step process involving several string intersection and annihilation events (hgure 
11). We expect that such a process probably will occur whenever the geometry is 
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Figure 11: Unlinking of two strings- the inverse of the process shown in hgure 
7- can occnr in several steps if the string tensions pnll in the right direction to nnlink 
the strings. A linking followed by two annihilations has the net resnlt of removing 
the short intervening segment and nnlinking the two longer strings. In this hgure, 
the basepoint is assumed to he behind the page, so that the dehnition of an s-string’s 
hux changes when it passes in front of any f-string. 
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appropriate for the unlinking of two strings, so that it is not necessary to perform the 
unlinking “by hand” in a single step within the simulation. 

(iv) Rearrangement: If two nodes approach each other closely but are prevented 
from annihilation by flux conservation reqnirements, several outcomes are conceivable. 
They may stick together and form a stable junction of more than three strings. They 
may bounce back and move apart again under the influence of string tensions. It is 
also plausible that the nodes could rearrange their connections and nndergo a sort of 
qnasi-intercommutation, as shown in hgure 12. In this simulation, we allow the nodes 
to bounce by introducing a small repulsion at distances shorter than dann- We also 
allow rearrangement with some probability, and examine the consequences of setting 
that probability either to zero or to some non-zero valne. The resnlts presented in 
section 6 indicate that rearrangement speeds up the network’s decay. 

We do not inclnde nodes of more than three strings as fnndamental objects, bnt 
it is qnite possible for a pair to become stnck together very close to each other. 
The string tension and short distance repulsion allow them to oscillate at a short 
distance, and snch a conhguration can behave like a single junction of fonr strings. 
Such adhesion becomes important under certain conditions, as we will see in section 
6 . 

4 Initial Conditions 

In order to perform onr dynamical simnlation, we must start with an initial 
conhguration. The generation of initial conditions models the symmetry-breaking 
transition which produces the strings. As is freqnently done, we use a lattice Monte 
Carlo procedure to generate an initial string network. The lattice spacing is to be 
identihed with the correlation length of the Higgs held at the time of string formation. 
The Higgs VEV is thus uncorrelated over distances longer than a lattice spacing, and 
at each site of the lattice, it takes a random value within the vacunm manifold. With 
a snitable interpolation along the length of each link, any plaquette of the lattice is 
mapped to some closed loop on the vacuum manifold. If this path is one of the non- 
contractable loops, then a string must pierce the plaquette. Each link of the plaquette 
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Figure 12: Rearranging the connections of two nodes which cannot annihilate. 
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is associated with an element of G which relates the Higgs values at the two ends of 
the link. The product of these elements around a closed plaquette must he within 
the unbroken group H, and can be taken as the flux of the string which pierces that 
plaquette. Strings which pierce the faces of a given unit cube must be joined together 
inside the cube in some way. If only two faces of the cube have non-trivial flux, then 
we interpret this as a single string segment passing through the cube. If three faces 
are non-trivial, we infer that there is a single vertex inside the cube. Cubes pierced by 
more than three ends require a more complicated arrangement of nodes and strings 
inside the cube. There may be more than one consistent way to join the string ends, 
and one must be chosen arbitrarily. 

Two different lattice Monte Carlo algorithms have been used for the current 
simulation to generate different initial distributions of strings. The first, very simple 
way to generate a random network of strings is to use an inhnite temperature lattice 
gauge theory: simply assign a random element of the unbroken group H = to each 
link of the lattice, and evaluate the product of links on the plaquette to hnd the flux 
through the plaquette. We refer to this as the lattice gauge Monte Carlo. There is 
no direct reference to a Higgs held in this technique. 

The other method we use is a discrete Higgs simulation analogous to that of 
references [27,28]. The essence of this method is that a discretized vacuum manifold 
is used. The breaking of continuous group G is modelled by using some discrete 
subgroup Q <Z G which contains the unbroken group H. Each lattice site is then 
assigned randomly to one of the discrete cosets in the space Q/H, corresponding to 
a choice of vacuum. With each link of the lattice there is associated an element 
of Q which transforms the Higgs field value at one end of the link to the value at 
the other. The element relating one coset to another is not unique; the possible 
elements themselves form a coset. The convention in this discrete Higgs method is 
to choose the “smallest” possible element for each link variable. “Small” is dehned 
with reference to a metric on the continuous group G\ if all elements are written in 
the form g = exp(iQ;T) where T is a normalized element of the Lie algebra of G, then 
the smallest element is the one with the smallest number a. In this way the Higgs 
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field is effectively interpolated in the smoothest possible way between lattice points. 

A suitable gauge transformation can be performed so that all Higgs field values 
he in the same coset, and all link variables lie within H, allowing all subsequent 
computations to be performed in terms of only H link variables. This is a lattice 
implementation of unitary gauge. 

For the present simulation, we take Q to be one of the discrete subgroups of SU (3). 
The simplest choice is the 24-element “dihedral-like” subgroup of SU{3) known as 
A(24)!^^' This group is generated by the matrices 



1 

0 

0 


(0 1 ox 


0 0-1 


0 0 1 


^ 0 -1 0 ) 


0 

0 


and 

A( 1 , 0 ) = diag(l, - 1 , - 1 ). 

A(24) is the smallest subgroup of SU (3) that contains S^, and in fact it is isomorphic 
to 5 * 4 , the permutation group on 4 elements: an isomorphism may be dehned which 
maps 

Ti ^ = (23), S+^s+ = (123) 

these generate the subgroup Hq ~ S '3 of permutations on elements 1-3. Hq may be 
viewed as the little group of a “Higgs” vev which has the form ho = (0, 0, 0,1): per¬ 
mutations of the hrst three positions leave ho invariant. The elements A(1,0) = 
diag(l, — 1 , — 1 ), A( 0 , 1 ) = diag(— 1 , — 1 , 1 ) and A(l, 1 ) = diag(— 1 , 1 , — 1 ) may be 
mapped to 5*4 by 

A(l, 0) ^ (14)(23), A(0,1) ^ (24)(13), A(l, 1) ^ (34)(12). 

These act nontrivially on ho, and together with the identity they generate the four 
distinct left cosets Ho, A{l,0)Ho, A{0,l)Ho, and A{l,l)Ho. Each of the three non¬ 
trivial cosets consists of the set of elements which transform ho to one of three other 
possible vev’s. For example, elements of A(l,0)iJo take (0,0,0,1) to (1,0,0,0). 
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To generate a network of strings, the discretized Higgs vev is randomly assigned to 
one of its four values at each lattice site. For a smooth interpolation, the link between 
two neighboring sites is chosen to be the smallest within the appropriate coset. In the 
coset H(l, 0)Ho, for example, the two smallest elements are H(l, 0)T2 (14) (23) (13) 

and H(1,0)T3 (14)(23)(12). Both of these have equal measure; a random choice 

may be made between them. In the identity coset, of course, the identity element 
is the smallest. Each other coset has two smallest elements of the form AT, where 
A G {H(1,0),H(0,1),H(1,1)} and T is one of the two transpositions in S'3 that fail 
to commute with A. 

After assigning vacuua and group elements, one can then transform to the unitary 
gauge in which the Higgs vev is the same at each site, and all flux information is 
encoded in Ss variables on the links, just as it is in the lattice gauge Monte Carlo. 

Properties of the initial network 

Both of the Monte Carlo algorithms described above create inhnite branched 
networks. The lengths of string segments between branching are distributed ex¬ 
ponentially, reflecting a constant probability of branching per unit length. In this 
respect, the two methods are similar, but the resulting networks differ in other sta¬ 
tistical properties. In table 1, we summarize some of these features. For comparison, 
we also include the corresponding information for the Z3 system (including both a 
Z3 lattice-gauge method and the tetrahedral discrete-Higgs simulation of [28]. The 
Z3 discrete-Higgs numbers are from reference [28] and the Z3 lattice gauge data are 
based on the author’s simulations (see also [30]). For each method, the fraction of 
plaquettes pierced by strings of each type is reported. (In the Z3 case, of course, 
there is only one type.) Note that in the lattice gauge method, each group element 
is weighted equally; therefore 2/3 of all plaquettes are pierced by strings in the Z3 
lattice gauge case, and 5/6 for Ss- Below this are the fractions of cubic lattice cells 
with 0,2,3,4,5 and 6 of their faces pierced by strings, the number densities of the 
different types of nodes per unit volume, the average length between junctions and 
the branching length (obtained from the exponential decay of the length distribution) 
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for each string species.* 

Pictures of typical initial string networks (figures 13) illustrate qualitatively the 
comparison between the different S '3 initial conditions. It is evident that for both Z 3 
and S '3 systems, the lattice gauge method produces a denser network: more cube faces 
are pierced by strings and more cells have high numbers of strings emerging through 
their faces. Correspondingly, in the lattice gauge method, fewer strings continue 
through more than one lattice cell without branching. 

Regardless of the initial conditions, stt junctions outnumber sss junctions in the 
S '3 network, and we can infer that a majority of the s strings end on a t string at at 
least one end. 

Although we have not extracted detailed statistics on the presence of disconnected 
loops and hnite networks in the initial distributions, they do not appear to form a 
signihcant part of the system. In this respect the S '3 networks are like the ones. 

The results presented in the next section show that in the S '3 model, unlike the 
Z 3 , differences in the statistical properties of the initial networks can have a very 
pronounced effect on the network’s evolution. 

5 Dynamical Results 

In this section, we present results from the dynamical simulation. The aim is 
to give a qualitative picture of the evolution and to determine which factors are 
most important in determining the fate of the network. The qualitative nature of 
the evolution depends on the ratio of the different string tensions. Two contrasting 
regimes are of interest: one in which the even strings are light, and another in which 
they are heavy. We also examine the influences of other factors. The results show a 
surprisingly marked dependence on statistical properties of the initial network. 

* Initially, string lengths are naturally clustered near integer multiples of the lattice spacing. 
The distribution looks smooth and exponential only when string lengths are placed in bins of 
at least one lattice spacing. For this reason, the decay length is not necessarily identical to the 
average segment length. 
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Method 

Zg L.G. 

Zs Higgs 

Ag L.G. 

^3 Higgs 

Faces with 



mmm 

■■1 

s string 

0.67 

0.52 



t string 

- 

- 

mm 

mm 

Cubes with 





0 ends 

0.01 

0.04 

0.00 

0.06 

2 ends 

0.12 

0.34 

0.01 

0.32 

3 ends 

0.17 

0.20 

0.05 

0.21 

4 ends 

0.38 

0.32 

0.21 

0.31 

5 ends 

0.24 

0.09 

0.20 

0.02 

6 ends 

0.09 

0.02 

0.55 

0.08 

Density of 
sss nodes 

0.56 

0.28 

0.22 

■ 

stt nodes 

- 

- 

1.36 

■HI 

Av. length 
bet. junctions 





s string 

2.38 

3.71 

1.02 

1.27 

t string 

- 

- 


1.59 

Branching length 



■H 


s string 

1.65 

3.33 



t string 

- 

- 

■a 

1.53 


Table 1: Statistics of Initial Conditions. 
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Lattice Gauge 



S String 
T String 


Figure 13: Sample Monte Carlo initial configurations (actually shown after a 

single time-step of dynamical evolution). The even, or 3-cycle, strings are shown in 

the darker color. Volume shown is 8^ in lattice units. 
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This section is organized as follows: First, we review some features of the behavior 
of Abelian, Z 3 strings, for comparison with the S '3 results. The Z 3 results are taken 
in part from [14] and in part from simulations by the author. The self-similar scaling 
evolution and the associated power-law behavior are demonstrated. After this come 
the S '3 results. We describe the dependence of these results on parameters of the 
simulation, focusing especially on two constrasting regimes of string tensions, and 
then discussing other important factors. Finally, we make a few comments about 
fluctuations and the effect of the hnite simulation volume. 

The Z 3 Network and Power-Law Evolution 

Since S '3 contains Z 3 as a subgroup, our program can easily be modified to simulate 
a Z 3 network by generating only s-strings in the initial conditions. Figure 14 shows 
the total string length and total number of nodes as a function of time for a typical 
run on a 30^ simulation volume. The initial conditions were generated with a Z 3 
lattice gauge method, rather than the discrete Higgs method of [14]. (In fact, these 
data were obtained in a simulation without any intercommutations — strings were 
allowed to pass through each other. There is no topological obstruction to prevent 
this in an Abelian network, and the inclusion of intercommutations makes only a 
small difference in the results.) The units used for this plot are the ones in which the 
initial lattice spacing is 1 and the string tension is T = 1. 

The time variable plotted on the x-axis in the figure is iAt, where i is the number 
of elapsed time steps. (This will be our convention for all remaining plots.) All 
distances and lengths are measured in units of the original lattice spacing. 

A transformation of the data shows more clearly the “scaling” behavior of the 
network. The scaling hypothesis is that the gross properties of the network are 
described by a single length scale, i, which grows with time as the network relaxes. If 
i is the typical distance between nodes, then the number of nodes per unit volume is 
n ~ In a scaling solution, the average length of a string segment between nodes 
is also ~ £, while the number of segments per unit volume is proportional to the node 
density, Hence the total string length per unit volume (or the energy stored in 
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strings) scales as l~‘^. In figure 15, we plot two scale variables with the dimensions 
of length; The inverse cube root of the number of nodes per unit volume, which we 
denote and the inverse square root of the string length per unit volume, Di. 
From hgure 28, we can see that both of these length scales are approximately equal 
and grow linearly with time. The slope of approximately 0.3 is close to that observed 
in reference [14]. Although the results plotted in £g. 28 and those of [14] were 
obtained from different initial conditions, the results agree very closely. Evidently, all 
noticeable differences disappear after just a few time steps. As with Z strings, the 
late-time evolution is essentially independent of the initial conditions. 

Evolution of the Z 3 network is self-similar in the sense that network looks sta¬ 
tistically the same at all times except for the increase in scale. The distribution of 
string segment lengths, for example, is exponential at all times, with only the scale 
changing. This can be seen in the semilogarithmic histograms of figure 27. Qualita¬ 
tively, a portrait of part of the network at a given time is indistinguishable from a 
suitably magnified portrait of an earlier time. 

5*3 network: general comments 

In the remainder of this section we describe results of the 5*3 simulation, which 
was run using a variety of different combinations of initial conditions, string tensions, 
and other parameter choices. Networks generated by lattice-gauge initial conditions 
were run on an 8 ^ simulation volume, while the less dense discrete-Higgs-generated 
networks were evolved on a volume of 12 ^. Statistical variables were obtained from 
averages over several runs. Computation time constraints made it unfeasible to run 
many times on larger volumes, but a few runs were performed on both larger and 
smaller volumes in order to examine the effects of finite size. 

Key results are displayed graphically in figures 17 - 27 and table 2 . Some results 
for an Abelian, Z 3 network are also shown for comparison. Figures 17 - 21 show three- 
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dimensional pictures of simulated cosmic string networks during their evolution. 
The pictures show several of the different patterns of evolution that can occur under 
different conditions. Pictures of a Z% network are provided as well. The plots in 
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D(n) 


Node Scale 




Figure 28: Scaling behavior of Z 3 strings on a 30^ volume. The typical distance 
between nodes, D^, grows linearly with time, and so does the inverse square root of 
the string energy density, Di. 
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figures 22 - 24 show the evolution of some length-scale variables as functions of time. 
These variables include the inverse cube root of the density of nodes, the average 
segment length between junctions for each type of string, and the two quantities Dg 
and Df. The last two are analogous to Di dehned above; the inverse square root 
of the string length per nnit volume, computed separately for each type of string. 
(Some of these variables may be rednndant.) With some exceptions discussed below, 
these length scale variables exhibit (after some transients at early times) the linear 
increase characteristic of self-similar evolntion. Such plots were made for simulations 
run under a wide variety of conditions giving a survey of the simulation’s parameter 
space. A few representative plots are shown here in order to show their typical shapes. 
The remainder are snmmarized in Table 2, which gives their slopes. The last set of 
hgures in this section, hgures 25 - 27, consists of a series of histograms which show how 
the distribution of string segment lengths evolves with time. Distribntions are shown 
for two contrasting cases discussed below, and corresponding data is also provided 
for the Z 3 network for the sake of comparison. 

The parameter space surveyed in Table 2 includes two different initial condition 
simulations, several different ratios of the string tensions, and two other binary choices 
affecting string collisions and close enconnters between nodes. The choice between 
“bridge” and “zipper” conhgurations for colliding non-Abelian strings is one choice 
(abbreviated B and Z in Table 2.) The other choice determines what happens when 
two neighboring nodes are within distance dann but are topologically unable to anni¬ 
hilate. In the cases labelled “R,” snch nodes nndergo a rearrangement of connections 
(see previons section) with a probability of 0.2 per time step. In the cases labelled 
“N,” no rearrangements are allowed, and the nodes simply bonnce. 

As a rough measure of the importance of intercommutations and NCC’s, Table 
2 also gives, in the last two columns, the ratio of the total number of NCC events 
to the number of intercommutations, and the ratio of the nnmber of NCC’s to the 
cnmnlative net nnmber of nodes annihilated. (For the Z 3 network, there are no NCC’s 
and the number of intercommutations is shown instead.) 


40 



Unless otherwise specified, all numerical quantities are given in units such that 
the initial lattice spacing is 1, the displacement of each node during a time step is 
given by Ax = At, and the magnitude of the tension Tf of the odd strings is 

normalized to 1. (This normalization was chosen because the special status of odd 
strings: they are the ones which cannot end, and in none of the cases simulated did 
they show a tendency to become extinct. Even strings, on the other hand, disappear 
almost completely under certain conditions.) The tension Tg of the even strings was 
varied relative to this fixed value. 

One salient feature of the results is hardly unexpected: the ratio of the tensions 
of different string species has a strong effect on the behavior of the network. The case 
of heavy s strings, in particular, is an exceptional one which will be discussed further 
below. Much more surprising is the strong difference in behavior between networks 
with different initial conditions. Evidently, different initial conditions lead to quite 
different trajectories which appear self-similar. 

We will address the issues in the following order: First, we discuss the different 
evolution patterns that occur with different choices of string tension, focusing on the 
contrasting limits of heavy and light even strings. Then we comment on other effects, 
including the effects of different initial conditions. Finally, at the very end of the 
section, we will make a few remarks about uncertainties and finite-size effects. 

S's network with light even strings 

Consider a network in which the even, or 3-cycle, strings have a much lower 
tension than the odd strings. In this case, the odd strings may pass through each other 
with comparative ease by forming new segments of the lighter even string. The odd 
strings may shrink with comparatively little energy cost in the creation of new strings. 
Furthermore, a zipper-type collision of two odd strings may have a result which, from 
the point of view of the odd string subsystem, resembles an intercommutation (see 
figure 16). The newly formed light string offers comparatively little resistance to the 
straightening of the rejoined heavy strings. We might expect that in this limit, the 
odd strings might behave as a network of Z 2 -like (i.e., unbranched and unoriented) 
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strings moving through a viscous medium formed by the branched Z 3 network of even 
strings. The odd strings can shorten and cut themselves into decaying loops while 
transferring part of their energy to the even string network. The even network, in 
turn, can dissipate its energy in much the same way as a Z 3 network. 

Simulation results seem consistent with the above picture. Figure 17 shows a 
series of portraits of the evolution of a network with tensions Tg = 0.5, Tt = 1.0, 
beginning from lattice gauge initial conditions. The conhguration at time t = 0 is in 
fact the one shown shown in £g. 13. In the early stages of evolution, we see a large 
number of highly crumpled, apparently Brownian t strings with thick webs of s string 
stretched between them. As time progresses, the odd strings begin to straighten and 
some closed loops shrink away and vanish. However, even though the odd strings are 
shortening at the expense of stretching the even ones, the population of even strings 
is at the same time being reduced through the annihilation of nodes, with the result 
that neither species of string becomes entirely dominant over the other. Instead, 
the evolution appears self-similar, with an approximately constant ratio between the 
amounts of the two different string species. This constant ratio is evident in figure 
23 from the linear evolution of both length scale variables Dg and Dt- Apparently, 
after some transient behavior at very early times, a sort of dynamical equilibrium is 
established, with energy being transferred in a steady cascade from the t network to 
s network and then lost to damping. 

Interestingly, decreasing the even string tension still further to 0.25 does not allow 
the odd strings to contract more quickly (see Table 2) - evidently they are impeded 
by the larger population of even strings. 

5*3 network with heavy even strings 

A situation which contrasts with the case of light s strings is one in which the s 
strings have a large tension. In this subsection, we discuss simulation results for the 
case with Tg = 2 , = 1 . A glance at figures 18 and 22 reveals that this leads to quite 

different results. The s, or even, strings are in this case only marginally stable against 
decay into pairs of t strings. If a heavy s string ends on a t string, then the two t 
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Figure 16: Two colliding t strings form an s-string zipper. If the s string is 
ignored, this event looks very much like an intercommutation. 
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segments pulling against the s string at the junction cannot prevent the s segment 
from shrinking unless the angle between the t segments is sufficiently small; in this 
marginal case of = 2 they must be collinear in order for the tensions to balance. 
In either lattice gauge or discrete Higgs initial networks, the majority of s strings end 
on t strings, and so there is very little to prevent these heavy strings from shrinking 
rapidly. 

The shrinking of an individual s segment frequently brings together a pair of 
nodes which are topologically unable to annihilate. The nodes stick together because 
to separate them again would require stretching the s segment. Similar adhesion 
happens whenever two non-commuting t-strings collide. Soon the network consists 
predominantly of odd strings stretched between small tangled clusters, and no fur¬ 
ther annihilation can occur until separate clusters meet, or the geometry changes 
sufficiently to allow some strings to pull free. Figure 22 shows the scaling variables as 
functions of time. The rapid disappearance of s string is evident from the extremely 
fast increase in the variable Dg in the upper plot, and from the lower plot which 
shows that the average length of an s segment quickly drops to near zero. This is an 
exceptional case: in all other cases besides Tg = 2, the average segment length for 
both string types increases with time. 

Even though the s strings do not obey the usual standard scaling behavior, it 
appears that the average t segment length does increase linearly with time, as do Dn 
and Dt- Evidently, the overall system obeys scaling once most of the s strings have 
shrunk. 

The prevalence of clusters connected by short s segments is apparent in the snap¬ 
shots of hgure 18 and is illustrated in a different manner by hgure 25, which shows 
semilogarithmic histograms of the distribution of segment lengths at a sequence of 
times. The length distributions for the two string types evolve in different ways. For 
lengths larger than one lattice spacing, the t segments follow an exponential distribu¬ 
tion (which appears linear on semilogarithmic axes) with a steadily increasing decay 
length. While the few long s segments that remain are also exponentially distributed. 
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there is no clear tendency for the decay length of this distribution to increase with 
time, and the number of long segments quickly shrinks into insignihgance compared 
with the large population of very short segments in the lowest bin of the histogram. 
The presence of so many short segments is associated with the adhesion of many pairs 
of nodes. 

Compare £g. 25 with the length distributions in the Tg = 0.5 case, shown in 
figure 26. In the latter, we see that both distributions develop a peak at short lengths 
(indicating that some clusters do occur), but the peaks are less pronounced and the 
number of such short segments decreases at a rate commensurate with the remaining 
population of strings. The exponential decay lengths for both distributions increase 
with time, as expected in a self-similar evolution. Evidently, the pile-up of clusters is 
less significant in the light-s-string case. Length distributions for the Abelian network 
(hgure 27) show no excess at all at short distances. 

Other Effects 

The non-Abelian simulation has been run under a variety of conditions. For the 
most part, the results are consistent with some type of self similar evolution. Variables 
having the dimensions of length all increase linearly with time. The exception is for 
networks with heavy s strings, in which the average s segment length decreases. This 
is the only case for which one type of string dominates completely over the other. In 
all other cases, the amounts of the two string types approach a constant ratio. 

The evolution of the network has some interesting dependence on parameters 
other than the string tensions. 

Most surprising is the very strong dependence on the statistics of the initial 
conditions. When the Higgs initial conditions are used, the network decays much 
faster, especially in the case of heavy s strings. Not only the absolute sizes, but 
the ratios of coefficients can have quite different values depending on which initial 
conditions are used. This dependence is somewhat mysterious. The lattice gauge 
initial conditions are much more dense than the Higgs, but the overall density of 
nodes decreases during the evolution, soon reaching a density comparable to that of 
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the initial Higgs-generated network, yet the two continue evolving at very different 
rates. Evidently, some statistical property other than the total density is important, 
such as, for example, the distribution of voids, or even correlations of the fluxes. A 
better understanding of this phenomenon should be the goal of future work. 

Notice that the details of the behavior of closely approaching nodes makes some 
difference in the evolution. The inability of some pairs of nodes to annihilate is 
apparently an important impediment to the network’s destruction. Rearrangements 
(quasi-intercommutations) of the connections of neighboring nodes evidently increase 
the mobility of flux, and increase the likelihood that eventually some neighboring 
pairs of nodes will be able to annihilate. As is apparent from Table 2, the inclusion 
of rearrangements almost invariably speeds up the decay of the network, often by 
roughly a factor of two. The choice of either zipper or bridge configurations for 
colliding non-commuting strings generally makes a smaller difference, if any. 

Finite-Size Effects, Fluctuations and Uncertainties 

A few remarks are in order concerning the interpretation of results and the effects 
of simulating on a finite volume. The size of our simulation volume is limited in 
practice by computation time; the flux computations (see Appendix) are computa¬ 
tionally intensive. Typical plots of scale variables like those in figures 22 - 24 exhibit 
some transient behavior at early times, followed by a period of linear increase. As 
the typical distance between nodes increases, however, a point is eventually reached 
where only a few nodes remain within the simulation volume. At this point, both 
systematic and random deviations from linear scale growth begin to occur. With 
fewer nodes, curves become bumpy due to lower statistics. Many strings wind all the 
way around cycles of the periodic boundary, so that the simulated network can no 
longer be expected to mimic an infinite one. A typical run meets one of two distinct 
fates at the end of its period of self-similar evolution. It may undergo a very sudden 
collapse as the last remaining nodes annihilate, leaving either no strings at all or a 
few strings stretched across the entire volume. This sudden collapse is visible as a 
sharp upward turn in the plots of length scale versus time. Alternatively, the network 
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0.25 

0.40 
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0.91 

Ts = 2.0, H 







■H 

B N 

2.7 

0.55 

0.69 

< 0 

1.2 

4.2 


Z N 

2.8 

0.46 

0.59 

< 0 

1.1 

4.5 

■il 

B R 

3.0 

0.71 

0.91 

< 0 

2.0 

2.9 

0.17 

Z R 

3.7 
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0.23 
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0.30 
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0.06 
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0.09 
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0.54 
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0.61 
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0.04 
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0.26 
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Z N 

■n 
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0.47 

0.43 
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0.20 

2.5 

0.06 
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0.44 

0.39 

0.29 

0.16 

0.17 
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0.05 
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tsi 

0.09 

0.17 

0.08 

0.04 

0.06 

1.9 

0.13 

Ts = 1, ^3 








(with intercommutations) 

0.37 

- 

0.38 

0.25 

- 

- 

0.12 


Table 2: Coefficients for network’s evolution: rates of change for Dg,Dt,Dn, and 
average segment length of each type, ratio of NCC’s to intercommutation events, and 
cumulative number of NCC’s per node, reported for different conditions. The hnal 
row gives Z 3 results. Abbreviations: LG=lattice gauge initial conditions, H=Higgs 
initial conditions, B=bridge NCC, Z=zipper NCC, R=rearrangements allowed, N=no 
rearrangements, IC=intercommutation, NCC=non-commutative collision. 
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- S Strings 

Strings 


Figure 17a: A series of snapshots showing the evolntion of an S '3 network with 
light strings. The heavier, odd strings are shown in darker color. We can observe 
long strings straightening and loops of t string contracting. The initial network was 
generated by the lattice gange Monte Carlo algorithm: it is the one shown in hgure 
13. The fnll 8 ^ simnlation volnme is shown. In the hrst frame, at t = 2, the odd 
strings are qnite crnmpled and are connected by a dense web of even strings. 
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Figure 17b: Light s strings, t = A 
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Figure 17c: Light s strings, t = Q 


50 









Figure 17d: Light s strings, t = 8 
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Strings 
- S Strings 


Figure 18a: Snapshots from the evolution of an network with heavy S strings, 
from lattice gauge initial conditions. The S strings, or even strings, are shown in 
thicker lines. The full 10^ simulation volume is shown. Notice the rapid shrinking 
of even segments, which leads to the formation of clusters that are slow to untangle. 
The hrst frame shown here is at f = 1.5. 
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Figure 18b: Heavy s strings, t = 3.0 
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Figure 18c: Heavy s strings, t = 4.5 
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Figure 18d: Heavy s strings, t = 6.0 
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Strings 
- T Strings 


Figure 19a: Evolution of a network from Higgs initial conditions with heavy S 
strings. The qualitative behavior is rather similar to that of the denser lattice gauge 
network, except that the network disappears much more rapidly, with large voids 
opening up very quickly. This hgure: t = 1.0. 
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Figure 19b: Heavy s strings, discrete Higgs initial conditions, t = 2.0 
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Figure 19c: Heavy s strings, discrete Higgs initial conditions, t = 3.0 
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S Strings 
Strings 


Figure 20a: Evolution of a network with equal string tensions, from lattice gauge 
initial conditions. First figure: t = 2.0. 


59 
















Figure 20b: Equal tensions, t = 4.0 
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Figure 20c: Equal tensions, t = 6.0 
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Figure 20d: Equal tensions, t = 8.0 
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Figure 21a: Snapshots from the evolution of a pure Z 3 network, included for 
comparison with the non-Abelian network. The hrst picture, shown here, is at t = 0. 
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Figure 21b: Z 3 network, t = 2.0 
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Figure 21b: Z 3 network, t = 4.0 
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Figure 21b: Z 3 network, t = 6.0 
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Figure 22: Scaling-law variables as a function of time for S '3 network with Tg = 2, 
lattice gauge initial conditions, bridge NCC and no rearrangement. 
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Figure 23: Length scales for S '3 network with light S strings. Tg = 0.5, lattice 
gauge initial conditions, bridge NCC, no rearrangements. 
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Figure 24: Length scales for S '3 network with equal string tensions, Higgs initial 
conditions, bridge NCC, rearrangement allowed. 
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Log N 


odd segment^ t=0.6 


Log N 


even segment^ t=0.6 




Figure 25: Evolution of the distribution of length of even (S') and odd (T) string 
segments between junctions, for a network with heavy S strings. 
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Log N 


t=1.2 


Log N 


t=2.4 



Figure 27: Segment length distributions for the Abelian Z 3 network, included 
here for comparison with hgures 25and 26. 
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may reach a stable or metastable configuration with a small number of nodes (typi¬ 
cally of order ten or fewer.) The characteristic length scale ceases its linear increase 
and reaches a plateau. Needless to say, we cannot make reliable inferences about 
an infinite network once the scale of our simulation volume becomes the important 
one, and we can only extract information on scaling behavior from the linear part of 
the curve. The conclusion that our data are consistent with scaling is based on the 
existence of a linear stretch which typically continues until at least one of the length 
scale variables has grown to one half the width of the simulation volume. We cannot 
rule out slow (e.g., logarithmic) deviations from scaling. 

With respect to fluctuations, not all simulation runs are alike. Runs with lattice 
gauge initial conditions and Tg = 2 seemed to undergo quite a bumpy evolution char¬ 
acterized by periods of temporary freezing alternating with cascades of annihilation. 
The bumpiness was smoothed out only by averaging over multiple runs. Runs with 
Tg = 2 without rearrangement almost always end with a few nodes remaining rather 
than with complete collapse. This is reflected in figure 22, where the curves of Dt 
and Dn versus time begin levelling off as the average segment length approaches 4 
(which is half the width of the volume). This levelling is probably a finite-size effect 
and not indicative of the behavior of an infinite network for the following reasons: It 
occurs when only a few nodes are left in the simulation volume, and it occurs at later 
times as the simulation volume is increased. In fact, one run performed on a larger 
10^ volume did not level off at all, but collapsed entirely instead. Higgs- generated 
networks tend to evolve more smoothly and to be more likely to undergo total collapse 
at the end of the simulation. 

6 Conclusions and Discussion 

One of the motivations for undertaking a simulation of non-Abelian string dynam¬ 
ics was to test for deviations from the familiar power-law behavior of the network’s 
energy density as a function of time, and especially to look for evidence for or against 
the conjecture that the tangling of strings would cause a non-Abelian network to freeze 
into a static (fixed in comoving coordinates) equilibrium. The results we have found 
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suggest that such a scenario, if it can occur at all, requires very special conditions. 
Over a range of different regimes, the results found here are consistent with some form 
of self-similar evolution, and with a density decaying at a rate commensurate with 
that of matter, n = The coefficient C, however, can vary in qnite interesting 

ways dne to non-Abelian effects. In a cosmological scenario, this coefficient controls 
the fractional contribution of strings to the energy of the universe. 

The interplay between the different string species causes the self-similar evolu¬ 
tion to be realized in some novel ways which can be qnite different nnder different 
conditions. Featnres such as the ratio of the populations of different string species 
may depend in rather complicated ways on various factors. Particularly striking and 
new is the strong inflnence of the initial distribntion of strings on the snbseqnent 
evolntion. Evidently, different self-similar evolntionary trajectories are possible, and 
initial conditions may be attracted to one trajectory or another depending on some 
statistical featnre other than the network’s overall density. This behavior seems al¬ 
most paradoxical, going against the notion that a scaling evolntion is one which has 
no memory of its initial state. 

There are interesting qnestions to investigate at both the “microscopic” and 
“macroscopic” level. By microscopic questions, we mean those concerning the dy¬ 
namics of individnal strings and individual collisions. Onr results indicate that the 
behavior of nodes in the network as they enconnter each other has a controlling in¬ 
fluence on the network’s evolution. The inability of non-Abelian nodes to annihilate 
is an important impediment to the removal of strings, and conhgnrations with nodes 
close together or coalescing into composites appear to be especially important in the 
interesting case of heavy s strings. The process we have called “rearrangement” of two 
nearby nodes seems to increase the mobility of strings and allow the network to decay 
more qnickly. This simulation was rnn with rather ad hoc assnmptions abont node 
collisions. A better microscopic understanding of nodes’ behavior and that of multi¬ 
node tangles will help provide inpnt for improved simnlations and nnderstanding of 
macroscopic questions. 
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At the macroscopic level, it would be desirable to obtain a better understanding 
of the principles governing the interplay between different string types and the strong 
influence of initial conditions. The types of network behavior seen here may open up 
new possibilities for cosmological model-building, aside from the string domination 
scenario. Some of the types of structures seen here, such as webs of light string 
stretched between heavier ones (as in the light-s-string scenario of this simulation) or 
the tangled clusters seen the Tg = 2 simulation might have interesting cosmological 
effects. We have seen hints that the strength of fluctuations in string density is 
different under different circumstances, being especially strong for Tg = 2. Might 
these fluctuations be interesting sources of structure formation? Along with analytical 
study, more rehned simulations of branched and non-Abelian networks may well prove 
worthwhile. 


APPENDIX A 

Simulation of non-Abelian Fluxes 

The subtle nature of non-Abelian magnetic flux in discrete gauge theories has 
been the subject of quite interesting research|^*’^'*’^°’^^’^^^ and presented algorithmic and 
computational challenges for this simulation. Ambiguities in the dehnitions of non- 
Abelian fluxes, mentioned in section 2, require that we design a very careful procedure 
to £x consistent dehnitions and maintain their consistency as the strings move. This 
must be done in order to allow us to make the appropriate comparisons of huxes when 
two strings collide or two nodes attempt to annihilate. Additional subtleties occur 
as a result of our using periodic boundary conditions: we are dealing with a discrete 
gauge theory on a non-simply-connected space manifold, and holonomies associated 
with non-trivial cycles become important. 

This appendix describes our procedure for dehning and comparing the huxes of 
non-Abelian strings. Further details may be found in [18]. As in section 2, we use 
the formalism developed in [16, 20]. 
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Gauge Fixing Conventions 

In our algorithm, the strings and nodes exist inside a rectangular volume with 
opposite sides identihed: a 3-torus. The subtleties associated with the periodic bound¬ 
ary conditions will be discussed later: for now we simply consider a network inside a 
rectangular volume with boundaries. We choose a cubic volume with one corner at 
(0, 0, 0) and with side length L. As explained in section 2, it is necessary to dehne 
fluxes using paths that begin and end at some basepoint. We choose a basepoint 
xq at the center of our simulation volume, xq = {L/2, L/2, L/2). Let each node be 
associated with a straight line segment (a “tail”) along the direction BN from the 
basepoint to the node’s location. Then let the flux of each outgoing string be de¬ 
fined with respect to a path which runs outward along this tail to a point which is 
taken to be vanishingly close to the the node. The path then encircles the string in 
a counterclockwise direction and returns to the basepoint along the node’s tail. This 
is illustrated in figure 29. This will be our convention for dehning the fluxes of the 
strings which join at a given junction. 

As illustrated in hgure 30, flux conservation requires that the product of all three 
fluxes emanating from a node be trivial when the fluxes are multiplied in a clockwise 
order with respect to the direction BN. i.e., if the strings in clockwise order are a, b, 
and c, then 

cba = e. (Al) 

In our algorithm, a record is maintained of the geometry of each node: the strings 
carry labels indicating the appropriate clockwise orientation. 

In the case of a doubly linked pair of nodes (hgure 5) two segments are collinear, 
and the order is therefore ill-dehned. In such a case we allow the order to be arbitrary, 
but the huxes of the two strings must be dehned in a way consistent with that order, 
such that the product of all three huxes is as usual trivial. The ordering must also 
be compatible between the two nodes which the segments join, so that the hux of a 
given segment is consistent at its two ends. (The consistency of segments from one 
end to another will be discussed below.) 
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Figure 29: Conventions for measuring the fluxes of the three strings emanating 
from a node. Each string’s flux is dehned as the flux through a path which leaves 
the basepoint xq along a straight line toward the node, then encircles the string in a 
counterclockwise direction as seen from the far end of the string and returns to the 
basepoint. 
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Figure 30: The composition 'f[3a of all three paths can be continuously deformed 
to a point. Therefore cba, the product of all three fluxes taken in a clockwise direction 
as seen from above the node, must be trivial. 


78 



The collection of standard paths dehned above represents a set of generators for 
7 ri(Af — {-D}). The flux state of a network of strings is fully specihed when we know 
the fluxes enclosed by all of these standard paths. The condition (Al) supplies one set 
of relations among these generators. For each string segment, there is also a relation 
involving the fluxes dehned at its two endpoints, as discussed below. 

“Sliding” hux from the endpoint 

By the conventions above, the hux of each string is dehned at its two endpoints. 
But for the purposes of this simulation it is necessary to make comparisons of the 
huxes of strings at arbitrary points along their lengths. For example, if two strings 
cross each other, their huxes must be compared at the crossing point in order to 
determine whether they commute. A meaningful comparison of the huxes of nearby 
string segments can be obtained only if the paths used to dehne those two huxes 
remain close to each other everywhere except in the immediate vicinity of the strings 
to be compared. In particular, the “tails” of the paths must not pass on opposite 
sides of any string, because such paths would give diherent hux measurements for the 
same string. It is possible to dehne the hux of a string at an arbitrary point along its 
length by sliding the standard path to the one which encircles the string at the point 
we wish to measure, as illustrated in hgure 31. If another string with hux b pierces 
the triangle which is swept out by the sliding path, then the hux at the new position 
is conjugated by b. If multiple strings occur, then the new hux a' is given by faf~^, 
where /, the total hux inside the triangle, is dehned as the product of the huxes of 
all enclosed strings, taken in order of increasing angle from the initial ray BP. The 
hux of each other string at the point where it pierces the triangle must in turn be 
dehned by a similar sliding procedure from one of its ends. This procedure, applied 
recursively, can thus dehne the hux of any string at an arbitrary point P along its 
length, as measured by a path which follows a straight line from xq towards P and 
encircles the string near P. 

If one slides the path all the way to the far end of the string, the resulting value 
of the hux must be consistent with the value measured by the standard path at 
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the other end. This specihes an additional set of relations among the generators of 
7 ri(A^ — {T*}) and furnishes one way of testing for errors in the simulation, as well 
as being necessary in order to dehne the fluxes of strings at the newly created nodes 
that result from string collisions. 

A modihed version of this sliding procedure is used to dehne all huxes initially 
from the original lattice. Paths composed of lattice links are deformed by a series of 
steps to straight-line paths from the base point to the location of each node. 

Holonomy Interactions 

As the network evolves dynamically and nodes change their position, the huxes 
dehned by these conventions may change in several diherent ways. First, as a node 
moves, its tail may be dragged across another string segment. Conversely, a string 
segment may be dragged across the node’s tail by the motion of other nodes. In both 
cases, the huxes of all strings at the node must be conjugated by the hux which is 
crossed, as shown in hgure 32. In addition, the geometry of the strings at a given 
junction may change, resulting in holonomy interactions among the three strings 
joined at that node. Such a process is shown in hgure 33: the motion of string a 
causes its standard hux to change, and also changes the clockwise ordering of the 
strings a, 6 ,and c. This requires both an adjustment of the hux dehnitions and of the 
order labels. 

Periodic Boundary Conditions Due to the periodic boundary conditions, our simu¬ 
lation volume has the topology of a 3-torus, T 3 . The non-trivial topology introduces 
three additional classes of noncontractible closed loops other than the ones associated 
with strings: namely those which wrap around one of the boundaries. These loops 
may be associated with nontrivial hux. As in the case of vortices on a Riemann 
surfacel^"^' the huxes of these loops may have holonomy interactions with the huxes 
of strings, and therefore a full description of string huxes requires us to maintain a 
record of these large loops. 

As representatives of the three “wraparound” classes, we choose canonical 
straight-line paths parallel to the coordinate axes, which we will refer to as P^;, P^, and 
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Figure 31: When the flux a of a string has been defined according to a path 
which encircles it near one end, the flux a' of that string at another point along its 
length can be defined by “sliding” the standard path a to a' as shown. If no other 
strings pierce the triangle which is swept out, then this merely represents a continuous 
deformation of a, and thus a' = a (fig. 31A). However, if the triangle is pierced by 
string with flux b as measured by path /?, then the flux is conjugated by b: a' = bab~^ 
(31B). More generally, if the triangle (or the oriented path 0 shown in 31C) encloses 
flux /, then a is conjugated by the total flux /, i.e., a' = faf~^. The total flux is 
given by the product of individual string fluxes, taken in order of increasing angle 
from the initial tail. (This can be seen by deforming a product of loops to a single 
loop enclosing all strings.) 81 





Figure 32: Crossing of a node’s tail by a string. This can happen either when 
the moving node drags its tail across the string (A), or when the string is dragged 
across the tail dne to the motion of another node (B). In both cases, the fluxes of all 
strings attached to the node whose tail is crossed must be conjugated by the flux of 
the crossing string. 
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Figure 33: Example of a holonomy interaction between strings attached to the 
same node. When the string carrying flux a is lifted over the other string carrying 
flux b, its flux must be redefined as a' = bab~^, and the conventional clockwise order 
of the three strings changes, with a and b exchanging places. The flux conservation 
condition is maintained: if cba = e originally, then also ca'b = e. 
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Figure 34: Transformation from one description of a flux to another at the 
boundary. Here a string is shown intersecting the plane of the page precisely where it 
intersects the boundary of the cubic simulation volume (dotted line). Under periodic 
boundary conditions, the two points labeled xq are identihed. The flux of the string 
may be described in terms of a path whose tail extends to the right of xq (a) or to 
the left {a'). If no other strings are present, then a is homotopically equivalent to 
In the more general situation shown in (B), a ^a'(0Lr0i?), and 

so the two descriptions of the flux are related through conjugation by /lC fR, where 
C is the flux associated with the path T and fi and /r are the overall fluxes enclosed 
by (pL and pR, respectively. The latter can be dehned in terms of paths lying entirely 
on one side or the other of the boundary. 
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r' 


(A) 



(B) 


Figure 35: Interaction between a string and one of the large loops of the 3-torus. 
As a string with flux a as dehned by the path a crosses to the right of xq as shown in 
(A), the path ah, where F is a straight-line path which wraps around the 3-torus, can 
be continuously deformed to the new straight-line path F'. Thus the flux C associated 
with F must be multiplied from the right by the string’s flux, C' = Ca. If the string 
crosses to the left as in (B), aV is deformed to F' and so the multiplication is from 
the other side: C = aC. 
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r^. Fj;, for instance, leaves the basepoint along the +£ direction, wraps aronnd the 
bonndary from a; = L to a; = 0, and then returns to the basepoint from the —x side. 
Along with values of the fluxes for all strings, our algorithm maintains a record of 
the As holonomies Cx,Cy, and Cz associated with F^,, Fy, and F^ respectively. These 
values must be known in order to make consistent comparisons of string fluxes across 
the boundary. (The procedure for doing so is illustrated by hgure 34.) The values 
of the Ci may change if a string crosses one of the curves F* (see figure 35). 
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